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ABSTRACT 



We report on the serendipitous discovery in the Blanco Cosmology Survey 
(BCS) imaging data of a z = 0.9057 galaxy that is being strongly lensed by a 
massive galaxy cluster at a redshift of z = 0.3838. The lens (BCS J2352-5452) 
was discovered while examining i- and z-band images being acquired in October 
2006 during a BCS observing run. Follow-up spectroscopic observations with 
the GMOS instrument on the Gemini South 8m telescope confirmed the lensing 
nature of this system. Using weak plus strong lensing, velocity dispersion, cluster 
richness A^ocb an d fitting to an NFW cluster mass density profile, we have made 
three independent estimates of the mass M200 which are all very consistent with 
each other. The combination of the results from the three methods gives M 2 oo = 
(5.1 ±1.3) x 1O 14 M , which is fully consistent with the individual measurements. 
The final NFW concentration C200 from the combined fit is C200 = 5-4lYi- We have 
compared our measurements of M 2 oo and C200 with predictions for (a) clusters 
from ACDM simulations, (b) lensing selected clusters from simulations, and (c) 
a real sample of cluster lenses. We find that we are most compatible with the 
predictions for ACDM simulations for lensing clusters, and we see no evidence 
based on this one system for an increased concentration compared to ACDM. 
Finally, using the flux measured from the [Oil] 3727 line we have determined the 
star formation rate (SFR) of the source galaxy and find it to be rather modest 
given the assumed lens magnification. 

Subject headings: gravitational lensing: strong — gravitational lensing: weak - 
galaxies: high-redshift 



Introduction 



Strong gravitational lenses offer unique opportunities to study cosmology, dark mat- 
ter, galactic structure, and galaxy evolution. They also provide a sample of galaxies, 
namely the lenses themselves, that are selected based on total mass rather than luminos- 
ity or surface brightness. The majority of lenses discovered in the past decade were found 
through dedicated surveys using a variety of techniques. For example, the Sloan Digital 
Sky Survey (SPSS) data have been used to effectively select le ns candidates from rich clus- 
ters (IHennawi et al.l 120081 ) t hrough intermediate scale clusters (lAUam et al.ll2007t iLin et al. 
2009h to i ndividual galaxies faolto n et al.lfioosl: Iwillis et al.ll2006l) . Other searches using the 



CFHTLS flCabanac et al.ll2007D and COSMOS fields flFaure et al. 



20081 : IJackson et al.ll2008l ) 
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have yielded 40 and 70 lens candidates respectively. These searches cover the range of giant 
arcs with Einstein radii 9 bin > 10" all the way to small arcs produced by single lens galaxies 
with 9 bin < 3". 

In this paper we report on the serendipitous discovery of a strongly lensed z = 0.9057 
galaxy in the Blanco Cosmology Survey (BCS) imaging data. The lens is a rich cluster 
containing a prominent central brightest cluster galaxy (BCG) and has a redshift of z — 
0.3838. Cluster-scale lenses are particularly useful as they allow us to study the effects of 
strong lensing in the core of the cluster and weak lensing in the outer regions. Strong lensing 
provides constraints on the mass contained within the Einstein radius of the arcs whereas 
weak lensing provides information on the mass profiles in the outer reaches of the cluster. 
Combining the two measurements allows us to make tighter constra ints on the mass M 2 oo and 
the concentration C200, of an NFW (jNavarro. Frenk. fc Whitelll995l ) model of the cluster mass 
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)). In addition, 



if one has spectroscopic redshifts for the member galaxies one can determine the cluster 
velocity dispersion , assuming the clu ster is virialized, and hence obtain an independent 
estimate for M200 fjBecker et al.l 120071). Finally one can also derive an M^nn estimate from 
the maxBCG cluster richness iV^oo ({Hansen et al.ll2005l ; I Johnston et al.l 120071 ) . These three 
different methods, strong plus weak lensing, cluster velocity dispersion, and optical richness, 
provide independent estimates of M 2 oo (^200 is defined as the mass within a sphere of 
overdensity 200 times the critical density at the redshift z) and can then be combined to 
obtain improved constraints on M200 and C200 Measurements of the concentration from 
strong lensing clusters is of particular interest as recent publications suggest that they may 
be more concentrated than on e would expect from ACDM models (IBroadhurst fc Barkana 
20081 : lOeuri fc Blandfordll2009h . 



The paper is organized as follows. In § [2] we describe the Blanco Cosmology Survey 
Then in § [3] we discuss the initial discovery and the spectroscopic follow-up that led to 
confirmation of the system as a gravitational lens, the data reduction, the properties of the 
cluster, the extraction of the redshifts, and finally the measurement of the cluster velocity 
dispersion and estimate of the cluster mass. In §H]we summarize the strong lensing features 
of the system. In §[5] we describe the weak lensing measurements. In §[6] we present the results 
of combining of the strong and weak lensing results and the final mass constraints derived 
from combining the lensing results with the velocity dispersion and richness measurements. 
We describe the source galaxy star formation rate measurements in §[7] and finally in §[S]we 
conclude. We assume a flat cosmology with Qm = 0.3, £7a = 0.7, and H = 70 km s _1 Mpc" 1 , 
unless otherwise noted. 
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2. The BCS Survey 

The Blanco Cosmology Survey (BCS) is a 60-night NOAO imaging survey program 
(2005-2008), using the Mosaic-II camera on the Blanco 4m telescope at CTIO, that has 
uniformly imaged 75 deg 2 of the sky in t he SPSS qriz bands in pr eparation for cluster finding 



with the South Pole Telescope (SPT) fjVanderlinde et al.ll2010i ) and other millimeter- wave 



experiments. The depths in each band were chosen to allow the estimation of photometric 
redshifts for L > L* galaxies out to a redshift of z = 1 and to detect galaxies to 0.5L* at 
5a to these same redshifts. The survey was divided into two fields to allow efficient use of 
the allotted nights between October and December. Both fields lie near 5 = —55° which 
allows for overlap with the SPT. One field is centered near a = 23.5 hr and the other is at 
a = 5.5 hr. In addition to the large science fields, BCS also covers 7 small fields that overlap 
large spectroscopic surveys so that photometric redshifts (photo-z's) using BCS data can be 
trained and tested using a sample of over 5,000 galaxies. 



3. Discovery of the lens and spectroscopic follow-up 

The lens BCS J2351-5452 was discovered serendipitously while examining i- and z- 
band images being acquired in October 2006 during the yearly BCS observing run. The 
discoverer (EJB-G) decided to name it "The Elliot Arc" in honor of her then eight-year old 
nephew. Table [I] lists the observed images along with seeing conditions. Fig. [I] shows a gri 
color image of the source, lens and surrounding environment (the pixel scale is 0.268" per 
pixel). The source forms a purple ring-like structure of radius ~ 7.5" with multiple distinct 
bright regions. The lens is the BCG at the center of a large galaxy cluster. Photometric 
measurements estimated the redshift of the cluster at z ~ 0.4, using the expected g — r 
and r — i red sequence colors, and also provided a photo-z for the source of z ~ 0.7, as 
described below. We note that this cluster was first reported as SCS O J235138— 545253 in 



an independent analysis of the BCS data by iMenanteau et al.l (120101 ) where its remarkable 



lens was noted and they estimated a photometric redshift of z — 0.33 for the cluster. 

We obtained Gemini Multi-Obj ect Spectrograp h (GMOS) spectra of the source and a 



number of the neighboring galaxies (ILin et al.l 120071 ). We targeted the regions of the source 



labeled A1-A4 in Fig. [2], and photometric properties of these bright knots are summarized 
in Table [2j In addition we selected 51 more objects for a total of 55 spectra. The additional 
objects were selected using their colors in order to pick out likely cluster member galaxies. 
Fig. El shows the r — % versus % color-magnitude diagram (top plot) and the g — r vs. r — % 
color-color diagram (bottom plot) of the field. The blue squares in the bottom panel of Fig. [3] 
show the four targeted knots in the lensed arcs. The green curve is an Scd galaxy model 
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( IColeman. Wu. fc Weedmanlll980f ) with the green circles indicating a photometric redshift 
for the arc of z ~ 0.7. Note this is not a detailed photo-z fit, but is just a rough estimate 
meant to show that the arc is likely at a redshift higher than the cluster redshift. Highest 
target priority was given to the arc knots and to the BCG. Then cluster red sequence galaxy 
targets were selected using the simple color cuts 1.55 < g — r < 1.9 and 0.6 < r — % < 0.73 
(also shown in the bottom panel of Fig. [3]), which approximate the more detailed final 
cluster membership criteria described below in §3.21 Red sequence galaxies with i < 21.6 
(3"-diameter SExtractor aperture magnitudes) were selected, with higher priority given to 
brighter galaxies with i(3") < 21. Additional non-cluster targets lying outside the cluster 
color selection box were added at lowest priority. 

We used the GMOS R150 grating + the GG455 filter in order obtain spectra with about 
4600 - 9000 A wavelength coverage. This was designed to cover the [Oil] 3727 emission line 
expected at ~ 6300 A, given the photo-z estimate of ~ 0.7 for the circs £ls well as the Me; 
absorption features at ~ 7000 A (and the 4000 A break at ~ 5600 A) for the z ~ 0.4 cluster 
elliptical galaxies. 

We used 2 MOS masks in order to fully target these cluster galaxies (along with the 
arcs) for spectroscopy. Each mask had a 3600 second exposure time split into 4 900-second 
exposures for cosmic ray removal. We also took standard Cu-Ar lamp spectra for wavelength 
calibrations and standard star spectra for flux calibrations. All data were taken in queue 
observing mode. A summary of the observations is given in Table [U 



3.1. Data Reduction 



The BCS imaging data were processed using the Dark Energy Survey data management 



system (DESDM V3) which is under develop ment at UIUC/NCSA/Fermilab f lMohr et al. 



20081 ; iNgeow et al.l 120061 ; IZenteno et al.ll201ll ). The images are corrected for instrumental 
effects which include crosstalk correction, pupil ghost correction, overscan correction, trim- 
ming, bias subtraction, flat fielding and illumination correction. The images are then astro- 
metrically calibrated and remapped for later coaddition. For photometric data, a photomet- 
ric calibration is appli ed to the single-epoch a nd coadd o bject photom etry. The A str(9matic 
softw are^ SExtractor jBertin fc Arnoutslll996l ). SCAMP feertirJl2006 |) and SWarp ( jBertin et al 
20021 ) are used for cataloging, astrometric refinement and remapping for coaddition over each 
image. We have used the coadded images in the griz bands for this analysis. 



http: / /www. astromatic.net 
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The spectroscopic data were processed using the standard data reduction package pro- 
vided by Gemini that runs in the IRAF framework^. We used version 1.9.1. This produced 
flux- and wavelength-calibrated 1-D spectra for all the objects. Additional processing for 
the source spectra was done using the IRAF task apall. 



3.2. Cluster properties 



We adopt the procedure used by the maxBCG cluster finder (IKoester et al.ll2007al Jbl) to 
determine cluster membership and cluster richness and to derive a richness-based cluster mass 
estimate. We first measure N ga i, the number of cluster red sequence galaxies, within a radius 
1 h~ x Mpc (= 4.55') of the BCG tha t are also brighter than 0.4L* at the cluster redshift 
z = 0.38. From IKoester et al.l (j2007al ). 0.4L* corresponds to an i-band absolute magnitude 
M = —20.25 + 51og/i at z = 0, while at z = 0.38, 0.4L* corresponds to an apparent 
magnitude i = 20.5 (specific value provided by J. Annis & J. Kubo, private co mmunication) , 



after a ccounting for both K-correction and evolution (also as described in IKoester et al. 



2007af ). We apply this magnitude cut using the SExtractor i-band MAG_AUT0 magnitude, 
which provides a measure of a galaxy's total light. (Note the 3"-diameter aperture magnitude 
used earlier for target selection in general measures less light cf. MAG_AUT0, but is better 
suited for roughly approximating the light entering a GMOS slit.) We set the red sequence 
membership cuts to be g — r and r — % color both within 2a of their respective central values 
(g — r )o = 1-77 and (r — z)o = 0.65, where the latter are determined empirically based on 
the peaks of the color histograms of galaxies within 1 h^ 1 Mpc of the BCG. In applying 
the color cuts we use the colors defined by SExtractor 3"-diameter aperture magnitudes 
(this provides higher S/N colors compared to using MAG_AUT0), and for the uncertainty we 
define a = a/ a color + intrinsic! where <J co i or is the color measurement error derived from the 
SExtractor aperture magnitude errors, and cr,„ h .,„ sjr is the intrinsic red sequence color width, 
taken to be 0.05 for g — r and 0.06 for r — i (IKoester et al.ll2007al ). 



Carrying out the above magnit ude and color cuts, w e obtain an initial richness estimate 



N, 



gal 



0.156 h 



44. Then, as discussed in lHansen et al.l ( 20051 1 . we define another radius 

„90.l 



Mpc = 1.51 h" 1 Mpc (= 6.88' 
the BCG to obtain a final richness estimate ./V200 



and repeat the same cuts within rf 00 of 



55. Finally, using the weak lensing 
mass calibration of I Johnston et al.l ( 120071 ) for maxBCG clusters, we obtain a mass estimate 
M 200 = (8.794 x 10 13 ) x (A" 20 o/20) L28 /i" 1 M© = (4.6 ±2.1) x 10 14 (h = 0.7), where we 
have also adopted the fractional error of 0.45 derived by lRozo et al.l (120091 ) for this iV 2 oo-based 



2 http://www. gemini.edu/sciops/data-and-results/processing-software 



-7- 



estimate of M200 for maxBCG clusters. 



We note that iRozo et al.l ( 120101 ) apply a factor of 1.18 to correct the iJohnston et al. 



111 



(120071) cluster masses upwar d, in order to account for a photo-z bias effect that is detailed 



Mandelbaum et al.l (120081 ) . We have not applied this correction as it makes only a OAa 
difference, although we remark that the resulting mass M200 — 5.4 x 10 14 M does appear 
to improve the (already good) agreement with our other mass estimates below (see §3.41 and 

§5U>. 



Fig. [3] shows color-magnitude and color-color plots of all galaxies that have i < 21 
(SExtractor MAG_AUT0) and that are within a radius rfj = 1.51 h' 1 Mpc (= 6.88') of the 
BCG. Note we have extended the magnitude limit here down to i — 21, to match the 
effective magnitude limit of our spectroscopic redshift sample ( §3.31 below) In particular, we 
find 86 maxBCG cluster members for % < 21, compared to the earlier A^oo — 55 for i < 20.5 
(corresponding to 0.4L*). These member galaxies are shown using red symbols in Fig. [3] and 
their properties are given in Table [3j 



3.3. Redshift determinations 

The redshift extraction was carried out using the xcsao and emsao routines in the 



IRAF external package rvsao ( jKurtz fc Minklll998l ). We obtained spectra for the 55 objects 



that were targeted. Four of these spectra were of the source. Out of the remaining 51 
spectra we had sufficient signal-to-noise in 42 of them to determine a redshift. Thirty of the 
objects with redshifts between 0.377 and 0.393 constitute our spectroscopic sample of cluster 
galaxies. Fig. [5] shows the spatial distribution of galaxies within a 6' x 6' box centered on the 
BCG, with maxBCG cluster members, arc knots, and objects with spectroscopic redshifts 
indicated by different colors and symbols. Table |3] summarizes the properties of the 30 cluster 
member galaxies with redshifts, and Table H] summarizes the properties of the remaining 12 
spectroscopic non-member galaxies. In Fig. [5] we show four examples of the flux-calibrated 
cluster member spectra including the BCG. 

Examination of Table |3] and Table H] shows that our spectroscopic sample is effectively 
limited at i « 21, as 39 of the 42 non-arc redshifts have % < 21. Note that of the 30 
spectroscopically defined cluster members, 22 are also maxBCG members, while another 
7 lie close to the maxBCG color selection boundaries. Also, of the 12 spectroscopic non- 
members, none meets the maxBCG criteria except the faintest one (with i = 21.58). 

The redshift of the source was determined from a single emission line at 71 00 A which is 
present with varying signal-to-noise in each of the knots that were observed. We take this line 
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to be the [OII]3727A line which yields a redshift of 0.9057 ±0.0005. The four flux-calibrated 
source spectra are shown in Fig. |6j Knot A2 was observed under seeing conditions that were 
a factor of two worse than for the other three knots (see Tabled]). 



3.4. Velocity dispersion and cluster mass measurement 



We used the 30 cluster galaxies t o estimate the red shift and velocity dispersion of the 
cluster using the biweight estimators of iBeers et al.l ( Il990l ). We first use the biweight location 
estimator to determine the best estimate for cz. This yields a value of cz = 115151.1 ± 
241.1 km s' 1 which translates to a redshift of z c = 0.3838 ± 0.0008. We then use this 
estimate of the cluster redshift to determine the peculiar velocity v p for each cluster member 
relative to the cluster center of mass using 

(cz-cz c ) 

We determine the biweight estimate of scale for v p which is equal to the velocity dispersion 
of the cluster. We find a value for the velocity dispersion of cr c = 855^gg 8 km s _1 . The 
uncertainties are obtained by doing a jackknife resampling. The redshift distribution is 
shown in Fig. [71 The overlaid Gaussian has a mean of z c and a width of o c x (1 + z c ). The 
lines represent the individual peculiar velocities v p of the cluster members. 

We can use the estimated velocity d ispersion t o derive an estimat e for the cluster mass. 
We use the results of lEvrard et al.l (120081 ) (see also lBecker et al.ll2007l ) which relates M 2 oo to 
the dark matter velocity dispersion 

l/a 



200 



10 15 Mp 



(2) 

h{z) V CTl5 J 

where h(z) = H(z)/ 100 km s _1 Mpc _1 is the dimensionless Hubble parameter. The values 
of the parameters were found to be a 15 = 1082.9 ± 4 km s _1 and a = 0.3361 ± 0.0026 
(lEvrard et al.ll2008l ). Using the standard definition of velocity bias b v = CTgai/ffDM, where 
a gai is the galaxy cluster velocity dispersion, we can rewrite Equation |2] as 

bl /a M 2m = 10 15 M * ( ^V /Q , (3) 

h{z) \cr 15 J 

where the quantity bl/ a M 2 oo parameterizes our lack of knowledge about velocity bias. Sub- 
stituting in the measured values for a ga i we obtain b l J a M 2 oo = 5.79l^gg x 10 14 M Q . 



Bayliss et al.l (120 111 and references therein) discuss an "orientation bias" effect which 



causes an upward bias in the measured velocity dispersions of lensing-selected clusters, due 



-9 - 



to the higher likelihood of the alignment al ong the line of sight of the major axes of the 
cluster halos, which are in general triaxial. iBayliss et al.l (1201 1 ) estimate that on average 
this will result in the dynamical mass estimate being biased high by 19-20%, u sing the same 



relati on between M 2 qo and velocity dispersion as we have used (Eqn. [2] above; lEvrard et al. 



20081 ). Correcting for this orientation bias effect would result in bl M200 = 4.8 x 10 14 M C 



which is not a significant difference, as the change is well under la. We therefore do not 
apply this correction, but we do note that it would improve the already good agreement with 
our other mass estimates in §3.21 and §6.11 (assuming no velocity bias, b v = 1.) 



4. Strong Lensing Properties 

We use the coadded r-band image shown in Fig. [8] to study the strong lensing features 
of the system as it has the best seeing and hence shows the most de tail. To remove the 



contribution to the arc fluxes from nearby objects we used GALFIT (jPeng et al.l 120021 ) to 
model the profiles of these objects (galaxies and stars) and then subtracted the model from 
the image. This was done for all four bands griz. These subtracted images are used for all 
determinations of arc fluxes and positions. A number of individual knots can be observed 
in the system along with the more elongated features. For example it appears that knot Al 
is actually composed of two individual bright regions which are resolved by the Sextractor 
object extraction described below. Knot A2 also appears to have two components although 
these are not resolved by the object extraction so we treat them as one in the modeling. Even 
though the cluster is fairly massive we do not see evidence for additional arc-like features 
outside of the central circular feature. In this case we expect the mass of the lens to be well 
constrained by the image positions. 

We use the criteria that to obtain multiple images the average surface mass density 
within the tangential critical curve must equal the critical surface mass density E cr ^. The 
tangentially oriented arcs occur at approximately the tangential critical curves and so the 
radius of the circle 9 nrc traced by t he arcs provides a measurement of the Einstein radius 



9ein ( jNarayan fc Bartelmannl 119961 ) . The mass Mein enclosed with the Einstein radius is 
therfore given by 

M EIN = E crit 7r(A^iJv) 2 (4) 

Substituting for S crit gives 

MEIN = AG^° EIN (5) 
where D s is the angular diameter distance to the source, D\ the angular diameter distance 
to the lens, and D\ s the angular diameter distance between the lens and the source. These 
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values are D s = 1610 Mpc, D t = 1081 Mpc and D ts = 825 Mpc. 



To determine the Einstein radius we ran Sextractor (IBertin fc Arnoutslll996l ) on the r- 
band image. This identified eight distinct objects in the image. We used the coordinates of 
those eight objects and fit them to a circle. The radius of the circle gives us a measure of the 
Einstein radius. The Einstein radius we measure is Oein = 7.53 ± 0.25" which translates to 
39.5±1.3kpc. This yields a mass estimate of (1.5±0.1) xlO 13 M and a corresponding velocity 
dispersion (assuming an isothermal model for the mass distribution) of a = 694 ± 12 km s -1 . 

The magnification of the lens fi ens can be roughly estimated under the assumption that 
the 1/2-light radius of a so urce at redshi f t z ~ 0.9 is about 0.46" (derived from the mock 
galaxy catalog described in IJouvel et al.l ( 120091 )). The ratio of the area subtended by the 
ring to that subtended by the source is ~ 0.6 x (AR/5r), where R is the ring radius and 5r 
is the 1/2-light radius of the source. The 0.6 factor accounts for the fraction of the ring that 
actually contains images. This gives a magnification of fi ens = 39. 

To o btain a more quantitativ e value for the magnification we have used the PixeLens^] 
program (ISaha &: Williams! 12004] ) to model the lens. PixeLens is a parametric modeling 
program that reconstructs a pixelated mass map of the lens. It uses as input the coordinates 
of the extracted image positions and their parities along with the lens and source redshifts. 
It samples the solution space using a Markov Chain Monte Carlo method and generates an 
ensemble of mass models that reproduce the image positions. We used the Sextractor image 
positions ob tain ed above and ass i gned the parities according to the prescription given in 
Read! (120071 ) . In ISaha fc Williams! (12004! ) they note that if one uses pixels that are too large 
then the mass distribution is poorly resolved and not enough steep mass models are allowed. 
We have chosen a pixel size such that this should not be a problem. 



It is well known (see for example ISaha &: Williams! (12006! )) that changing the slope of 
the mass profile changes the overall magnification, in particular a steeper slope produces 
a smaller magnification but does not change the image positions. Therefore the quoted 
magnification should be taken as a representative example rather than a definitive answer. 
The magnification quoted is the sum over the average values of the magnification for each 
image position for 100 models. We obtain a value of fi ens = 141 ± 39 where the error is 
the quadrature sum of the RMS spreads of the individual image magnifications. PixeLens 
can also determine the enclosed mass within a given radius. For the 100 models we obtain 
M EIN = (1.4 ± 0.02) x 1O 13 M which is within la of the mass obtained from the circle fit 
described above. 



3 Version 2.17: 



http: / /www.qgd.uzh.ch/ programs/pixelens/ 
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In order to combine the strong lensing mass with the mass estimate from the weak 
lensing analysis (in §6. II below) we will need to estimate the mass within Oein that is due to 
dark matter alone (Mdm)- To do this we will need to subtract estimates of the stellar mass 
(Mg) and th e hot gas mass (Mq) from the total mass M EIN . To determine M s we use the 
GALAXEV (IBruzual fc Charlotl 120031 ) evolutionary stellar population synthesis code to fit 
galaxy spectral energy distribution models to the griz magnitudes of the BCG within the 
Einstein radius. The BCG photometric data are taken from the GALFIT modeling described 
above, and we sum up the light of the PSF-deconvolved GALFIT model inside the Einstein 
radius. The GALAXEV models considered are simple stellar population (SSP) models which 
have an initial, instantaneous burst of star formation; such models provide good fits to early- 
type galaxies, such as those in clust ers. In particula r we find a good fit to the BCG, using a 
SSP model with solar metallicity, a IChabrierl ( 120031 ) stellar initial mass function (IMF), and 
an age 9.25 Gyr (this age provided the best \ 2 over the range we considered, from 1 Gyr 
to 9.3 Gyr, the latter being the age of the universe for our cosmology at the cluster redshift 
z = 0.38). The resulting stellar mass (more precisely the total stellar mass integrated over 
the IMF) is M s = 1.7 x 10 12 M Q . 



To estimate the gas mass Mq we have looked at estimates of hot gas fraction /, 



gas 



i n cluster cores from X-ray observations. Typi cal f gas measurements are of order 10% 
( IMaughan et al.l 12004 iPointecouteau et al.l 120041 ) which give us an Mq estimate of 1.5 x 
10 12 M m . 



Finally we calculate the total M/L ratio within Oein for the i-band. This yields a value 
of (M/L)i = 33.7 ± 4.4 (M/L) Q . 



5. Weak Lensing Measurements 



5.1. Adaptive Moments 



We used the program Ellipto (ISmith et al. to compute adaptive moments (IBernstein fc Jarvisl 



20021 ; iHirata et al.ll2004f ) of an object's light distribution, i.e., moments optimized for signal- 
to-noise via weighting by an elliptical Gaussian function self-consistently matched to the 
object's size. Ellipto computes adaptive moments using an iterative method and runs off 
of an existing object catalog produced by SExtractor for the given image. Ellipto is also 
a forerunner of the adaptive moments measurement code used in the SDSS photometric 
processing pipeline Photo. 



We ran Ellipto on our coadded BCS images and corresponding SExtractor catalogs, 
doing so independently in each of the griz filters to obtain four separate catalogs of adaptive 



second moments: 
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where I(x,y) denotes the measured counts of an object at position x,y on the CCD image, 
and w(x,y) is the elliptical Gaussian weighting function determined by Ellipto. The images 
are oriented with the usual convention that North is up and East is to the left, i.e., right 
ascension increases along the —x direction and declination increases along the +y direction. 
We then computed the ellipticity components e\ and of each object using one of the 
standard definitions 



For each filter, we then identified a set of bright but unsaturated stars to use for PSF 
fitting. We chose the stars from the stellar locus on a plot of the size measure Q xx + Q yy 
from Ellipto vs. the magnitude MAG_AUT0 from SExtractor, using simple cuts on size and 
magnitude to define the set of PSF stars. We then derived fits of the ellipticities e±,e2 and 
the size Q xx + Q yy of the stars vs. CCD x and y position, using polynomial functions of cubic 
order in x and y (i.e., the highest order terms are x 3 , x 2 y, xy 2 , and y 3 ). On each image, these 
fits were done separately in each of 8 rectangular regions, defined by splitting the image area 
into 2 parts along the x direction and into 4 parts along the y direction, corresponding to the 
distribution of the 8 Mosaic-II CCDs over the image. This partitioning procedure was needed 
in order to account for discontinuities in the PSF ellipticity and/or size as we cross CCD 
boundaries in the Mosaic-II camera. Also note that the individual exposures comprising the 
final coadded image in each filter were only slightly dithered, so that the CCD boundaries 
were basically preserved in the coadd. To illustrate the PSF variation in our images, we 
present in Figure [9] "whisker plots" that show the spatial variation of the magnitude and 
orientation of the PSF ellipticity across our i- and r-band images . In addition, we also show 
the residuals in the PSF whiskers remaining after our fitting procedure, showing that the 
fits have done a good job of modeling the spatial variations of the PSF in our data. 




(9) 
(10) 



5.2. PSF Modeling 



13 



We next used our PSF model to correct our galaxy sizes and ellipticities for the effects 
of PSF convolution. Spec ifically, for the size measure Q xx + Q yy we used the simple relation 



(cf. Iffirata fc Seliakl 120031 ) 



Qxx,true + Qyy,true = {Q xx,observed + Qyy , observed) ~ (Q xx, PSF Qyy, PSF) (H) 

to estimate the true size Q XX) true + Q yy ,true of a galaxy from its observed size Q X x,observed + 
Qyy, observed) where Q X x,psf + Q yy ,PSF is obtained from the PSF model evaluated at the x, y 
position of the galaxy. For the ellipticities we similarly used the related expressions 

&i,observed ( 1 1 \ ■ 1 o /io\ 

ei,tme = d h 1 - -5- ei iPSjP , z = 1, 2 (12) 



-^2 V -^2 , 

D — 1 Qxx,PSF + Qyy, PSF / 1Q x 

2 = 7) IT) ^ ' 

^£ xx , observed 1~ ^£ yy , observed 

The relations used in this simple correction procedure strictly hold only for unweighted 
second moments, or for adaptive moments in the special case when both the galaxy and the 
PSF are Gaussians. We have therefor e also checked th e resul ts using the more sophisticated 



linear PSF correction" procedure of iHirata fc Seljakl (120031 ). which uses additional fourth 



order adaptive moment measurements (also provided here by Ellipto) in the PSF correction 
procedure. In particular, the linear PSF correction method is typically applied in weak 
lensing analyses of SDSS data. However, we found nearly indistinguishable tangential shear 
profiles from applying the two PSF correction methods, and we therefore adopted the simpler 
correction method for our final results. 



5.3. Shear Profiles and Mass Measurements 



Given the estimates of the true galaxy ellipticities from Equation f fl2|) . we then computed 
the tangential (e p ) and B-mode or cross (e x ) ellipticity components, in a local reference frame 
defined for each galaxy relative to the BCG: 



tT = t\ cos(20) — t2 sin(20) 
e x = ei sin(20) + e 2 cos(20) 



(14) 
(15) 



where <fi is the position angle (defined West of North) of a vector connecting the BCG to the 
galaxy in question. Here we have dropped the subscript true for brevity. The ellipticities 
were then converted to shears 7 using 7 = e/R, where R is the responsivity, for which we 
adopted the value R = 2(1 — a^ N ) = 1.73, using <tsn = 0.37 as t he intrinsic galaxy shap e 
noise as done in previous SDSS cluster weak lensing analyses (e.g., iKubo et al.ll2007l . 120091 ). 
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We then fit our galaxy shear measurements to an NFW profile by minimizing the fol- 
lowing expression for x 2 : 



E[lfi - iNFwjfh ^200, C 20 o)] S 
i=l 7 



(16) 

where the index i refers to each of the N galaxies in a given sample, r» is a galaxy's pro- 
jected physical radius from the BCG (at the redshift of the cluster), a 1 is the measured 
st andard deviation of the ga laxy shears, and ^nfw is the shear given by Equations (14-16) 
of I Wright fc Brainerdl (120001 ) for an NFW profile with mass M 2 oo and concentration 0200- We 



used a standard Levenberg-Marquardt nonlinear least-squares routine to minimize x 2 and 
obtain best-fitting values and errors for the parameters M200 and C200 of the NFW profile. 



ten been used to constrain the mass distributions of galaxy clusters (e.g., 


King & Schneider 


2001; 


Clowe & Schneider 


2001; 


Kubo et al. 


2009; 


Oeuri et al. 


2009; 


Okabe et al. 


2010h. Note 



that we chose the above expression for x 2 since it does not require us to do any binning in 
radius, but for presentation purposes below we will have to show binned radial shear profiles 
compared to the NFW shear profiles obtained from our binning-independent fitting method. 

For the shear fitting analysis, we defined galaxy samples separately in each of the four 
griz filters using cuts on the magnitude MAG_AUT0 and on the size Q X x,observed + Qyy,observed, 
as detailed in Table The bright magnitude cut was chosen to exclude brighter galaxies 
which would tend to lie in the foreground of the cluster and hence not be lensed, while the 
faint magnitude cuts were set to the photometric completeness limit in each filter, as defined 
by the turnover magnitude in the histogram of SExtractor MAG_AUT0 values. For the size 
cut, we set it so that only galaxies larger than about 1.5 times the PSF s ize would be used, 



as ha s been typically done in SDSS cluster weak lensing analyses (e.g., iKubo et al.l 120071 . 
20091 ). Note that in order to properly normalize the NFW shear profile to the measure- 
ments, we also need to calculate the critical surface mass density S crit , which depends on 
the redshifts of the lensed source galax i es as w ell as the redshift of the lensing cluster; see 
Equations (9,14) of IWright fc Brainerdl ( 120001 ). To do this, we did not use any individual 
redshift estimates for the source galaxies in our analysis, but instead we calculated an ef- 
fective value of 1/S cr j f via an integral over the source galaxy redshift distribution publi shed 



P 

for the Canada- France-Hawaii Telescope Legacy Survey (CFHTLS; lllbert et al.l 120061 ). as 
appropriate to the magnitude cuts we applied in each of the griz filters. 

Our NFW fitting results are shown in Figures ITDlfTTl and detailed in Table We show 
results for both the tangential and B-mode shear components. As lensing does not produce 
a B-mode shear signal, these results provide a check on systematic errors and should be 
consistent with zero in the absence of significant systematics. For all of our filters, our 
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B-mode shear results are indeed consistent with no detected mass, as the best-fit M 2 oo is 
within about la of zero. On the other hand, for the tangential shear results in the r, i, 
and z filters, we do indeed obtain detections of non-zero M 20 o at the better than 1.5<r level. 
In the g filter we do not detect a non-zero M 2 oo- Comparing the weak lensing results from 
the different filters serves as a useful check of the robustness of our lensing-based cluster 
mass measurement, in particular as the images in the different filters are subject to quite 
different PSF patterns, as shown earlier in Fig. |9j Though the mass errors are large, the M 20 o 
values from the r-, i-, and z-band weak lensing NFW fits are nonetheless consistent with 
each other and with the masses derived earlier from the velocity dispersion and maxBCG 
richness analyses. Moreover, independent of the NFW fits, we have also derived probabilities 
(of exceeding the observed \ 2 ) that our binned shear profiles are consistent with the null 
hypothesis of zero shear. As shown in Table El we see that the B-mode profiles are in all 
cases consistent with zero, as expected, but that the tangential profiles for the r and i filters 
are not consistent with the null hypothesis at about the 2a level (probabilities « 0.06), thus 
providing model-independent evidence for a weak lensing detection of the cluster mass. 



5.4. Combining Weak Lensing Constraints from Different Filters 

Here we will combine the weak lensing shear profile information from the different 
filters griz in order to improve the constraints on the NFW parameters, in particular on 
M 2 oo- The main complication here is that although the ellipticity measurement errors are 
independent among the different filters, the most important error for the shear measurement 
is the intrinsic galaxy shape noise, which is correlated among filters because a subset of the 
galaxies is common to two or more filters, and for these galaxies we expect their shapes to 
be fairly similar in the different filters. In particular we find that the covariance of the true 
galaxy ellipticities between filters is large, for example, the covariance of e\ between the % 
and r filters, Cov(ei 5 j, ei >r ) = XX e M ~~ 61,1) ( e i,r — &i,r), is about 0.9 times the variance of 
ei in the i and r filters individually. The same holds true for e 2 and for the other filters as 
well. We will not attempt to use a full covariance matrix approach to deal with the galaxy 
shape correlations when we combine the data from two or more filters. Instead, we take 
a simpler approach of scaling the measured standard deviation of the shear (the a 1 used 
to calculate x 2 m Equation [16]) by V 'N / N unique , where N is the total number of galaxies 
in a given multi-filter sample, and N un i que is the number of unique galaxies in the same 
sample. This is equivalent to rescaling \ 2 i n the NFW fit to correspond to N unique degrees of 
freedom instead of N. We have verified using least-squares fits to Monte Carlo simulations of 
NFW shear profiles that this simple approach gives the correct fit uncertainties on M 20 o and 
c 2 oo when the mock galaxy data contain duplicate galaxies, with identical e± and e 2 values, 
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simulating the case of completely correlated intrinsic galaxy shapes among filters. Note that 
our approach is conservative and will slightly overestimate the errors, because the galaxy 
shapes in the real data are about 90% correlated, not fully correlated, among filters. 

Before fitting the combined shear data from multiple filters, we make one additional 
multiplicative rescaling of the shear values, so that all filters will have the same effective value 
of corresponding to a fiducial effective source redshift z crit = 0.7. This correction is 

small, with the largest being a factor of 1.18 for the z-band data. The results of the NFW fits 
for the multi-filter samples are given in Table EJ where we have tried the filter combinations 
i + r, i + r + z, and i + r + z + g. We see that these multi-filter samples all provide better 
fractional errors on M 2 oo compared to those from the single-filter data. Also, as expected, the 
B-mode results in all cases are consistent with no detected M 20 o and zero shear. For our final 
weak lensing results, we adopt the NFW parameters from the i + r + z sample, as it provides 
the best fractional error (<tm 20 o/-^2oo ~ 0.5) on M 2 oo; we obtain M 20 o = 5.0^ 2 ; 3 x 1O 14 M , 
and c 2 oo = 4.9l 2 ' 2 . Figure [121 shows the shear profile data and best fit results for the i + r + z 
sample. This final weak lensing value for M 2 oo agrees well with the earlier values of M 2 oo 
derived from the cluster galaxy velocity dispersion (assuming no velocity bias) and from the 
cluster richness A^oo- 



6. Combined Constraints on Cluster Mass and Concentration 

6.1. Combining Strong and Weak Lensing 

In this section we combine the strong lensing and weak lensing information together in 
order to further improve our constraints on the NFW profile parameters, in particular on 
the concentration parameter c 2 oo- The addition of the strong lensing information provides 
constraints on the mass within the Einstein radius, close to the cluster center, thereby 
allowing us to better measure the centra l concentration of t he NFW profile and improve the 



uncertainties on the concentration c 2 oo- lOguri et al.l (120091 ) incorporated the strong lensing 



information in the form of a constraint on the Einstein radius due to just the dark matter 
distribution of the cluster, and they specifically excluded the contribution of (stellar) baryons 
to the Einstein radius. Their intent, as well as ours in this paper (§ I6.2p . is to compare the 
observed cluster NFW concentration to that predicted from dark-matter- only simulations. 
Thus the contribution of baryonic matter should be removed, most importantly in the central 
region within the Einstein radius, where baryonic effects are the largest due in particular to 
the presence of the BCG. In practice with the present data we can do this separation of the 
baryonic contribution only for the strong lensing constraint, and strictly speaking the weak 
lensing profile results from the total mass distribution rather than from dark matter alone. 
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Here we combine the strong a nd weak lens i ng da ta using an analogous but somewhat 



simpler method compared to that of lOguri et al.l ( 120091 ). specifically by adding a second term 



to x 2 (Equation [TB]) that describes the constraint on the dark matter (only) mass within the 
observed Einstein radius: 

2 _ hi ~ iNFwjrj] M 2 qo, c 200 )] 2 [M DM (< Oe) - M NFW (< E ] M 20 o, c 20 o)] 2 

i=l °7 U M DM (<0 E ) 

where Oe = 7.53" is the observed Einstein radius due to the total cluster mass distribution, 
Mdm(< Oe) is the dark matter (only) mass within 9e, and Mnfw(< Oe] M 200 , c 2 oo) is the 
mass within 9 E of an NFW profile with mass M 20 o, concentration c 20 o, redshift z = 0.38, and 
source redshift z = 0.905 7. M NF w(< F ] M 200 , c 200 ) is derived based on Equation (13) of 



Wright fc Brainerdl (120001 ). As obtained earlier in §U we estimate M DM (< E ) by subtracting 
estimates of the stellar mass and hot gas mass from the total mass within 9 E , obtaining 
Mdm(< 9e) = (1.18 ± 0.2) x 1O 13 M when subtracting off both stellar and gas mass, or 
M DM (< 9 E ) = (1.33 ± 0.2) x 1O 13 M when subtracting off only stellar mass. The former is 
our best estimate of M DM (< 0e), while the latter serves as an upper limit on M DM (< 9 E ) 
and hence on the best- fit concentration c 2 oo- We also conservatively estimate the error on 
Mdm(< 0e) to be one of the stellar mass/gas mass components added in quadrature to the 
uncertainty on the total Mein from £JU 

We apply the combined strong plus weak lensing analysis to our best weak lensing 
sample, the multi-filter i + r + z data set. The fit results are given in Table [5] and shown 
in Figure [12j We find M 2 oo = 4.9l 2 ' 2 x 10 14 solar masses, nearly identical to the final weak 
lensing result. We also get a concentration c 20 o = 5.5^ 2 ;g, again consistent with the final weak 
lensing fit, but with a 30% improvement in the error on c 2 oo, demonstrating the usefulness 
of adding the strong lensing information to constrain the NFW concentration. Using the 
upper limit Me>m{.< 0e) value (with only stellar mass subtracted) gives nearly the same 
M 2 qq = 4.8^2 2 x 1O 14 M , while the resulting NFW concentration is higher, as expected, 
with c 20 o = 6.2l 3 ; 2 , but still consistent with the fit using our best estimate of M DM (< 0e)- 



6.2. Combining Lensing, Velocity Dispersion and Richness Constraints 

In the above sections we have obtained quite consistent constraints on the cluster mass 
M 20 o using three independent techniques: (1) M 200 (lensing) = 4.9± 2 ; 2 x 1O 14 M from com- 
bined weak + strong lensing ffTTT): (2) M 200 (a c ) = 5.79tfH x 1O 14 M from the cluster 
galaxy velocity dispersion a c ( §3.4j assuming no velocity bias, b v = 1); and (3) M 20 o(iV 2 oo) = 
(4.6 ±2.1) x 1O 14 M from the maxBCG-defined cluster richness A^ 200 f %T2]) . We note that 
these methods are subject to different assumptions and systematic errors. For example, the 
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velocity dispersion based mass estimate assumes the cluster is virialized, an assumption sup- 
ported by the Gaussian-shaped velocity distribution of the cluster members shown in Fig. [71 



Also, the richness based mass estimate relies on the A^otr-^oo calibration ( iJohnston et al. 



20071 ) obtained for SDSS maxBCG clusters at lower redshifts z — 0.1 — 0.3 and assumes that 
this calibration remains valid for our cluster at z = 0.38. It is encouraging that we are ob- 
taining a cluster mass measurement that appears to be robust to these disparate assumptions 
and that shows good agreement among multiple independent methods. 

We will therefore combine the results from the different techniques in order to obtain 
final constraints on M200 and concentration C200 that are significantly improved over what 
any one technique permits. Specifically, we can add the M200 constraints from the velocity 
dispersion and richness measurements as additional terms to the weak + strong lensing x 2 
(Equation [T7|) : 

2 ^ hi ~ Twwfc; M200, c 2 qo)] 2 [M DM (< 6 E ) - M NFW (< d E \ M 2 qo, C200)] 2 

i=l °7 °Mdm(<0e) 

[M 200 (a c ) - M 200 ] 2 , [M 200 (N 200 ) -M200] 2 , 1Q , 



2 2 V / 

a M 200 (a c ) a M 20 o(N 2 oo) 

Minimizing this overall \ 2 results in the final best-fitting NFW parameters M200 = 5.lti3 x 
1O 14 M and C200 = 5.4l}' 4 . These results are consistent with the final lensing-based values 
M 2 oo(lensing) = 4.9^^2 x 1O 14 M and c 2 oo(lensing) = 5.5l 2 ;g, but have errors nearly a factor 
of two smaller. Note these quoted errors are 1-parameter, la uncertainties; we plot the joint 
2-parameter, la and 2a contours in Fig. \~ 



We also note that for the three methods weak lensing, velocity dispersion, and cluster 
richness, the corresponding NFW parameters result from the total mass distribution, con- 
sisting of both dark matter and baryonic (stellar plus hot gas) components. Dark matter 
is dominant over the bulk of the cl uster, while baryons can have a significant effect in the 



cluster core (e.g.. lOguri et al.ll2009f ). As described earlier (§ 16.11) . we have thus subtracted 
off the baryonic contribution to the strong lensing constraint as the intent is to compare 
(see below) our cluster concentration value against those from dark-matter-only simulations. 
Note that we have not isolated the dark matter contribution for the other three methods 
and cannot easily do so. For weak lensing, the shear profile is sensitive to the total mass 
distribution, not just to dark matter. For the velocity dispersion method, the galaxies act 
as test particles in the overall cluster potential, which is due, again, t o both dark matter 



and baryons. For the cluster richness method, the IJohnston et al.l ( 120071 ) A^oo-^oo relation 
we use was derived from stacked cluster weak lensing shear profile fits, including a BCG 
contribution but otherwise no other baryonic components; thus again the M200 value is es- 
sentially for the total mass distribution. Nonetheless, the bulk of the baryonic contribution 
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is in the cluster core and is accounted for via the strong lensing constraint, so we expect the 
comparison below of our cluster concentration value to those of dark matter simulations to 
be a reasonable exercise. 



Recent analyses (e.g., lOguri et al.l l2009t iBroadhurst fc Barkanal 120081 ) of strong lens- 
ing clusters have indicated that these clusters are more concentrated than would be ex- 
pected from ACDM predictions, th ough others have ar gued that no discrepa ncy exists if 



baryo nic effects are accounted for (IRichard et al.ll2010l ). In the former case, lOguri et al. 
(120091 1 found a concentration c w - lT ~ 9 for the 10 strong lensing clusters in their anal 
ysis sample, compared to a value of c w ~ 6 expected for strong-lensing-selected clus- 



ters or c vir « 4 for clusters overall (e.g., IBroadhurst fc Barkanal l2008t lOguri fc Blandford 



20091). We il lustra te these different concentration values in Fig. [T3J We use Eqn. (17) of 
Qguri et al.l (I2009I ). c vir (sim) = (l | 7 : 8 ) n 71 (M vir /2. 78 x lO 12 Af )-°- 081 , which comes from the 



ACDM N-body simulations of iDuffy et al.l (|2008|) . to show the t ypical concentration of clus- 
ters overall, and multiply by a factor of 1.5 (jOguri et al. 2009 ) to sho w the higher c oncen - 
tration expected for lensing selected clusters. We also use Eqn. (18) of lOguri et al.l ( 120091 ). 
c vir (fit) = (1 ^o.7i (Myir/lO^Me)" - 081 , to show the fit results for their cluster sample. In 

0.4 to match the redshift of our cluster. Moreover, we convert 



(i+^)"- 71 
these relations, we set z 

from the M vir , c vir convention used by Qguri et al. ( 2009]) to our M? n n , c 20 o convention, using 



the detailed rel a tions found in Appendix C of iHu fc Kravtsovl (120031 ) or in the Appendix of 
Johnston et al.l (120071 ). For the plotted M 20 o range, it turns out that C200 ~ 0.83 c vir . From 
Fig. H31 we see that our best-fit value of C200 = 5.4!}^ is most consistent with the nominal 
ACDM concentration value for lensing-selected clusters, and does not suggest the need for 
a concentration excess in this particular case. It's likely that larger strong lensing cluster 
samples will be needed to more robustly compare the distribution of concentration values 
with the predictions of ACDM models. 



7. Source Galaxy Star Formation Rate 



We can use the [OII]3727 line in the c alibrated spectra described in § 13.31 to estimate 
the star formation rate (SFR). As noted by iKennicuttl ( 19981 ) the luminosities of forbidden 
lines like [OII]3727 are not directly coupled to the ionizing luminosity and their excitation 
is also sensitive to abundance and the ionization state of the gas. However the excitation of 
[Oil] is well behaved enough that it can be calibrated through Ha as an SFR tracer. This 
indirect calibration is very useful for studies of distant galaxies because [OII]3727 can be 
observed out to redshifts z ~ 1.6 and it has been measured in several large samples of faint 
galaxies (see references in IKennicuttl (119981 ) ) . If we know the [Oil] luminosity then we can 
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use equation 3 from iKennicuttl (119981 ) to determine a star formation rate for the galaxy 

SFR(M yr" 1 ) = (1.4 ± 0.4) x 10- 41 (L[OII])(ergs s^ 1 ) (19) 

where the uncertainty reflects the range between blue emission-line galaxies (lower limit) 
and more luminous spiral and irregular galaxies (upper limit). 

As noted above, in order to extract the SFR we need to determine the total source flux 
from the [Oil] line. We determine this using 



[Oil] 



(20) 



where f{v)[oii\ is the total flux emitted by the source in the [Oil] line, f{y)h is the flux 
measured in the [Oil] line in each spectrum, f[y)s is the flux in the knot spectrum contained 
within the i-band filter band pass and j{y)i is the flux from the source in the z-band. 

Using the GALFIT-subtracted «-band image we determine j[y)i by summing the flux 
in an annulus of width 3" that encompasses the arcs. The flux j[y)h is measured by fitting a 
gaussian plus a continuum to the [Oil] line in each spectrum and integrating the flux under 
the gaussian fit. The flux f{v)s is calculated as follows. For each spectrum we first fit the 
continuum level, we then add the fitted continuum plus the [Oil] line flux and convolve it with 
the filter response curve for the SDSS z-band filter and integrate the convolved spectrum. 

We have determined f{y)\oii\ separately for each knot that was targeted for spectra. 
The fluxes are listed in Table |6]for each knot. We convert j{v)\on\ into an [Oil] luminosity 
and then use Equation [19] to determine a star formation rate for each knot. This rate is the 
raw rate which must be scaled by the lens magnification f\ ens to determine the true rate. 
We quote the SFR for t he two values of fi ens that were determined in We assume one 
magnitude of extinction (lKennicuttlll998l ) and have corrected the measured [Oil] luminosity 
to account for this. This yields the star formation rates listed in Table [6] for the two values 
of fiens- The rate for knot A3 is higher by a factor of 2 compared to the others because it 
has a small f(v)s compared to the other knots but the value of f{y)L is quite similar to 
the other knots. This can clearly be seen in Figure EJ We can combine the measurements 
for the four knots using a simple average to quote an overall SFR. This yields values of 
SFR(/ /ens = 49) = 4.6 ± 0.7 and SFR(/ iens = 141) = 1.3 ± 0.2. 



These rates are si gnificantly smal ler that those obtained for the 8 o'clock arc ( lAllam et al. 



20071 ) and the Clone (ILin et all 120091 ) which were 229M yr" 1 and 45M Q yr" 1 respectively 
(after converting to our chosen cosmology). Both these systems were at much higher red- 
shift (2.72 and 2.0 respectively) so one would potentially expect higher rates from these 
systems. They also had smaller values of fiens- We can compare ou r result to blue galax- 
ies at similar redshift from the DEEP2 survey fjCooper et al.ll2008l ). Using Figure 18 of 
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Cooper et al.l (120081 ) we obtain a median SFR of about 34M yr 1 for a redshift z = 0.9 



galax y which is also high er than our measurement. Other measurements using the AEGIS 
field (INoeske et al.l 120071 ) give a median SFR ranging from 10M Q yr~ x to 4OM yr -1 de- 
pending weakly on the galaxy mass, which is u nknown in ou r case. Our measurement can 
be compared to the far-right plot of Figure 1 in iNoeske et al.l ( 120071 ) and we fall on the low 
side of the measured data. Note that these conclusions are dependent on the magnification 
values used, for example smaller values such as those obtained for the Clone or the 8 o'clock 
arc would yield larger values for the SFR. 



8. Conclusions 

We have reported on the discovery of a star-forming galaxy at a redshift of z — 0.9057 
that is being strongly lensed by a massive galaxy cluster at a redshift of z — 0.3838. 

The Einstein radius determined from the lensing features is 0ein = 7.53 ± 0.25" and 
the enclosed mass is (1.5 ± 0.1) x 1O 13 M , with a corresponding SIS velocity dispersion of 
a = 694 ± 12 km s _1 . 

Using GMOS spectroscopic redshifts measured for 30 cluster member galaxies, we ob- 
tained a velocity dispersion a c = 855lgg 8 km s -1 for the lensing cluster. 

We have derived estimates of M 2 oo from measurements of (1) weak lensing, (2) weak + 
strong lensing, (3) velocity dispersion a c , and (4) cluster richness N 2 oo = 55. We obtained the 
following results for M 200 : (1) M 200 (weak lensing) = 5.0l 2 ;g x 1O 14 M , (2) M 20 o(lensing) = 
4.9+2:2" x 1O 14 M , (3) M 20 oO c ) = 5 - 79 ±L99 x 1O 14 M (assuming no velocity bias, b v = 1), and 
(4) M 2 oo(-^2oo) = (4.6 ±2.1) x 1O 14 M . These results are all very consistent with each other. 
The combination of the results from methods 2, 3 and 4 give M 200 = 5.ll};| x 1O 14 M , which 
is fully consistent with the individual measurements but with an error that is smaller by a 
factor of nearly two. The final NFW concentration from the combined fit is C200 = 5.4l}' l5 
which is also consistent with the lensing-based value but again with a smaller error. 

We have compared our measurements of M 2 oo and c 2 oo with predictions for (a) clusters 
from ACDM simulatio ns, (b) lensing^ selec ted clusters from simulations, and (c) a real sample 



of cluster lenses from lOguri et al.l (120091 ). We find that we are most compatible with the 
predictions from ACDM simulations for lensing clusters, and we see no evidence that an 
increased concentration is needed for this one sys tem. We are stud ying this further using 



other lensing clusters we observed from the SDSS (IDiehl et al.ll2009l ). These clusters will be 
the subject of a future paper. 
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Finally, we have estimated the star forming rate (SFR) to be between 1.3 to 4.6 M yr -1 , 
depending on magnification. These are small star-forming rates when compared to some of 
our previously reported systems, and are also small when compared with rates found for other 
galaxies at similar redshifts. However we caution that this conclusion is entirely dependent 
on the derived lens magnification. 
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Fig. 1. — A gri color image of the Elliot Arc and its cluster environment. The scale is 
indicated by the horizontal arrow. 
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Fig. 3. — (Top) r — i vs. i(MAG_AUT0) color-magnitude diagram for all galaxies (black points) 
with i < 21 and within a radius rf^ = 1.51 h~ l Mpc (= 6.88') of the BCG. Colors are 
measured using 3"-diameter aperture magnitudes. Galaxies meeting the maxBCG cluster 
color selection criteria (see §3.2() are plotted in red, with red circles indicating cluster members 
brighter than i = 20.5, and red squares indicating fainter cluster members. (Bottom) g—r vs. 
r — i color-color diagram for the same galaxies as in the top panel. Red circles and squares 
again indicate brighter and fainter maxBCG cluster members, while the black rectangle 
indicates the color selection box (approximating the more detailed maxBCG color criteria) 
used to select likely cluster galaxies for GMOS spectroscopy (see $3]). In addition, the 4 
bright knots A1-A4 (Fig. |2]) in the lensed arcs are shown by the blue squares. The green 
curve is an Scd galaxy model (IColeman. Wu. &: Weedmanl Il980l ) at redshifts z = — 2, 
with green circles highlighting the redshift range z = 0.65 — 0.75, indicating an approximate 
photometric redshift z ~ 0.7 for the arc knots. 
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Fig. 4. — Relative positions of all galaxies (points) with i(MAG_AUT0) < 21 within a 6' x 6' box 
centered on the BCG. Cluster member galaxies defined using maxBCG criteria (see §3.2p are 
plotted in red, with red circles indicating members brighter than i = 20.5, and red squares 
indicating fainter members. The 4 bright knots A1-A4 (Fig. |2]) in the lensed arcs are shown 
by the blue squares. Galaxies determined to be cluster members from GMOS redshifts are 
plotted with open magenta circles, while those found spectroscopically to be non-members 
are shown with open cyan triangles (see §3.31) . North is up and East is to the left. 
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Fig. 5. — Four examples of flux-calibrated cluster member spectra (in f u ). The spectra have 
been smoothed (with a boxcar of 5 pixels = 17.8 A) to improve the signal-to-noise ratio. 
The spectrum in the top left is that of the BCG. The prominent absorption features used in 
the redshift identification are marked. 
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Fig. 6. — Flux-calibrated spectra (in f u ) for the knots A1-A4. The spectra have been 
smoothed (with a boxcar of 5 pixels = 17.8A) to improve S/N. Knot A2 was observed under 
seeing conditions that were a factor of two worse than for the other three knots. The [O II] 
3727 A line is marked. 




Fig. 7. — The redshift distribution for the 30 cluster members in Table [3j The tick marks 
at the top represent the individual cluster member peculiar velocities. The solid line is a 
Gaussian with mean and sigma equal to z c and o c x (1 + z c ) respectively (see §3.4p . 




Fig. 8. — The coadded r-band image. The lensing features 
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Fig. 9. — (Top panels) "Whisker" plots that show the clear spatial variation of the PSF 
ellipticity vs. CCD x, y position in our i- (left) and r-band (right) images. The size of each 
whisker is proportional to the PSF ellipticity epsF = yjef PSF + e\ PSF , where a whisker with 
ellipticity e = 0.1 is shown at the top center of the figure. Each whisker is oriented at an 
angle 9p$F = \ taxi~ l (e2,psF / &i,psf) counterclockwise from horizontal. (Bottom panels) The 
corresponding whisker plots after subtraction of the PSF model described in §5.21 showing 
the removal of the bulk of the spatial variation of the PSF ellipticities. 
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Fig. 10. — The points with error bars show the tangential (top) and B-mode (bottom) radial 
shear profiles for the galaxy sample used for weak lensing analysis in the % (left) and r (right) 
filters. In each panel, the solid curve shows the shear profile for the best-fitting NFW mass 
density profile, as determined via the procedure described in §5.31 The dashed horizontal 
lines indicate zero shear. The best-fit NFW parameters and details of the galaxy sample are 
given in Table 
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Fig. 11.- 



Similar to Figure [TU1 but for the z (left) and g (right) filters. 
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Fig. 12. — Similar to Figure [IUJ but for the multi-filter i + r + z sample. For the tangential 
shear profile fits in the top panel, the long-dashed curve gives the results using weak lensing 
only, while the dotted and solid curves give the results using combined weak plus strong 
lensing. The dotted curve is for the case where we estimated the dark matter mass within 
the Einstein radius by subtracting off just a stellar mass contribution, while the solid curve 
is for the case where we also subtracted off an estimated gas mass contribution. See §5.4} 
§6.11 and Table |5] for details. 
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Fig. 13. — Confidence contours for the best-fitting NFW mass M 2 oo and concentration C2oo> 
obtained by combining the lensing, velocity dispersion, and cluster richness constraints, as 
described in §6.21 The 2-parameter, la contours are shown in solid blue, while the 2a 



contours are shown in hatched blue. The outer dashed contours show the 2-parameter, 2a 
constraints derived solely from the weak + strong lensing analysis of §6.11 Also, as described 
in §6.2} the 3 mostly horizontal curves show the concentration vs. mass relation at z = 0.4 
for: (bottom) clusters overall from ACDM simulations; (mid dle) lensing selecte d clusters 
from simulations; and (top) a real lensing cluster sample from lOguri et al.l (120091 ) . 
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Table 1. Observation Log 



BCS Imaging 

g 14 Dec 2006 2x125 sec 1.44" 

r 14 Dec 2006 2x300 sec 1.29" 

g 11 Nov 2008 2x125 sec 1.03" 

r 11 Nov 2008 2x300 sec 0.88" 

% 30 Oct 2006 3x450 sec 1.18" 

z 30 Oct 2006 3x450 sec 1.31" 

GMOS spectroscopy 

GG455 4 Aug 2007 4x900 sec 0.56" Mask 1 includes knots A1,A3,A4 

GG455 4 Aug 2007 4x900 sec 1.14" Mask 2 includes BCG and knot A2 

GG455 4 Aug 2007 1x5 sec - Cu-Ar Mask 1 

GG455 4 Aug 2007 1x5 sec - Cu-Ar Mask 2 

GG455 14 Aug 2007 1x5 sec - 1.5" slit 

GG455 14 Aug 2007 lx90sec 0.95" Standard star EG21 
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Table 2. Knots Targeted for Spectroscopy 



Knot RA a Dec a i(3") b g - r c r - i c 



Al 


357.912477 


-54.881691 


21.94 


0.85 


0.77 


A2 


357.911467 


-54.882801 


21.49 


0.81 


0.69 


A3 


357.906225 


-54.883464 


22.30 


0.84 


0.74 


Al 


357.907100 


-54.879967 


21.46 


0.91 


0.77 



a RA and Dec are epoch J2000.0 and are given in degrees. 

b i-band magnitudes for the knots are computed in 3"- 
diamctcr apertures, after first subtracting a model of the BCG 
light derived using the Galfit galaxy fitting program jPeng et al] 
l2002h . 

c g — r and r — i colors are computed from 3"-diameter SEx- 
tractor aperture magnitudes. 
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Table 3. Cluster Galaxies 



Object ID a 


RA b 


Dec b 


i(MAG_AUT0) a 




9 ~ r a 




r — i a 




rcdshift z c 


maxBCG Cluster Members' 1 


















15173 (BCG) 


357.908555 


-54.881611 


17.36 ±0.00 


1 


86 ± 0.01 





68 ± 


00 


0.3805 ± 0.0003 


16097 


357.972190 


-54.856522 


18.58 ±0.00 


1 


87 ±0.02 





71 ± 


.01 




16926 


358.069064 


-54.838013 


18.67 ± 0.01 


1 


.72 ± 0.02 





.75 ± 


01 




14954 


357.990606 


-54.881805 


18.70 ± 0.01 


1 


.78 ± 0.02 





.69 ± 


01 




14458 


357.922935 


-54.891348 


19.05 ± 0.00 


1 


.77 ±0.02 





.66 ± 


.01 


0.3844 ± 0.0002 


15111 


357.911305 


-54.879770 


19.10 ± 0.01 


1 


76 ± 0.03 





.64 ± 


01 




13772 


357.854114 


-54.908062 


19.21 ± 0.01 


1 


76 ± 0.02 





.65 ± 


01 




14873 


357.913389 


-54.883120 


19.22 ± 0.01 


1 


.78 ± 0.02 





.66 ± 


01 




15204 


357.917968 


-54.874795 


19.22 ± 0.01 


1 


.81 ± 0.02 





.64 ± 


01 


0.3827 ±0.0002 


15305 


357.929749 


-54.874524 


19.32 ± 0.01 


1 


.66 ± 0.02 





.66 ± 


.01 




15124 


357.915316 


-54.877257 


19.39 ± 0.01 


1 


.72 ± 0.03 





.65 ± 


.01 


0.3929 ± 0.0005 


11813 


357.856326 


-54.957697 


19.40 ± 0.01 


1 


.70 ± 0.03 





.63 ± 


01 




13629 


357.781583 


-54.911494 


19.57 ± 0.01 


1 


.67 ±0.03 





.64 ± 


.01 




16084 


357.932492 


-54.855191 


19.62 ± 0.01 


1 


.81 ± 0.03 





.67 ± 


.01 


0.3864 ± 0.0003 


14828 


357.858498 


-54.884039 


19.69 ± 0.01 


1 


.76 ± 0.03 





.67 ± 


01 




15056 


357.929263 


-54.878393 


19.69 ± 0.01 


1 


.83 ± 0.03 





.71 ± 


.01 




13028 


357.836914 


-54.923702 


19.69 ± 0.01 


1 


.70 ± 0.03 





.65 ± 


.01 




14267 


357.742239 


-54.897565 


19.70 ± 0.01 


1 


.66 ± 0.04 





.73 ± 


01 




13939 


357.743518 


-54.903441 


19.72 ± 0.01 


1 


78 ± 0.04 





70 ± 


01 




14892 


357.917045 


-54.881040 


19.75 ± 0.01 


1 


.67 ±0.03 





.58 ± 


.01 




17276 


358.061857 


-54.827190 


19.83 ± 0.01 


1 


.68 ± 0.03 





.67 ± 


01 




12997 


357.992988 


-54.925261 


19.85 ± 0.01 


1 


79 ± 0.04 





70 ± 


01 




14685 


357.948912 


-54.885316 


19.85 ± 0.01 


1 


.77 ±0.03 





.70 ± 


.01 




14727 


357.914364 


-54.884540 


19.86 ± 0.01 


1 


.74 ±0.04 





.64 ± 


01 




12907 


357.971817 


-54.926860 


19.88 ± 0.01 


1 


.73 ± 0.03 





.71 ± 


01 




15525 


357.891439 


-54.867148 


19.95 ± 0.01 


1 


,85 ± 0.04 





.67 ± 


.01 


0.3802 ± 0.0005 


13874 


357.767222 


-54.904760 


19.98 ± 0.01 


1 


.76 ± 0.04 





.68 ± 


.01 




14875 


357.896375 


-54.880676 


20.00 ± 0.01 


1 


.81 ± 0.04 





.67 ± 


01 




14827 


357.956282 


-54.883059 


20.02 ± 0.01 


1 


.73 ± 0.04 





.69 ± 


.01 




14169 


357.942454 


-54 896963 


20 05 + 01 


1 


69 + 04 


o 


.67 ± 


01 




14620 


357.906482 


-54.885446 


20.09 ±0.01 


1 


.65 ±0.03 





.66 ±0 


01 


0.3822 ± 0.0004 


11254 


357.792343 


-54.968633 


20.11 ±0.01 


1 


.76 ±0.03 





.63 ±0 


.01 




19279 


357.988291 


-54.784591 


20.12 ±0.01 


1 


.67 ±0.04 





.67 ±0 


.01 




15027 


357.880749 


-54.878200 


20.16 ±0.01 


1 


.76 ±0.04 





.63 ±0 


01 


0.3876 ± 0.0006 


12805 


357.947638 


-54.929596 


20.18 ±0.01 


1 


.70 ±0.04 





.76 ±0 


01 




13899 


358.003818 


-54.902294 


20.19 ±0.01 


1 


.84 ±0.04 





.66 ±0 


.01 




14741 


357.943783 


-54.884299 


20.21 ±0.01 


1 


.85 ±0.05 





.69 ±0 


01 




12671 


358.055413 


-54.931583 


20.22 ±0.01 


1 


.71 ±0.04 





.59 ±0 


01 




14843 


357.901141 


-54.880772 


20.23 ±0.01 


1 


.72 ±0.04 





.61 ±0 


.01 




14088 


357.936003 


-54.898050 


20.27 ±0.01 


1 


.79 ±0.05 





.69 ±0 


.01 


0.3816 ± 0.0005 


14969 


357.910388 


-54.878452 


20.29 ±0.01 


1 


.76 ±0.05 





.64 ±0 


01 




12875 


357.935888 


-54.927451 


20.30 ±0.01 


1 


.80 ± 0.04 





.64 ±0 


.01 




13537 


357.937414 


-54.911304 


20.31 ±0.01 


1 


.75 ±0.04 





.68 ±0 


.01 


0.3849 ± 0.0003 


15314 


357.902668 


-54.872019 


20.34 ±0.01 


1 


.70 ±0.05 





.60 ±0 


01 


0.3841 ± 0.0004 


14669 


357.916522 


-54.885626 


20.36 ±0.01 


1 


.89 ± 0.06 





.73 ± 


.02 


0.3862 ± 0.0003 
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Table 3 — Continued 





DAb 


uec 


•i(MAG_AUT0) a 




a 

g — r 


■a 
T — % 


rcQsliiit z 


14639 


357.954384 


-54.885672 


20.36 ±0 


01 


1 


.77 ±0.05 


0.67 ± 0.01 




14232 


357.904683 


-54.895183 


20.38 ±0 


01 


1 


.76 ±0.04 


0.64 ±0.01 


0.3882 ± 0.0003 


14703 


357.865075 


-54.884639 


20.41 ±0 


.01 


1 


,81 ±0.05 


0.62 ± 0.02 




14690 


357.914016 


-54.883857 


20.42 ±0 


.01 


1 


.68 ± 0.04 


0.64 ± 0.01 




15463 


357.882797 


-54.868348 


20.43 ±0 


01 


1 


.71 ±0.05 


0.63 ±0.01 


0.3877 ± 0.0004 


16005 


357.991736 


-54.856356 


20.44 ±0 


.01 


1 


,70 ±0.05 


0.66 ±0.01 




15333 


357.975175 


-54.872090 


20.44 ±0 


.01 


1 


,69 ±0.05 


0.63 ±0.01 




14972 


357.909534 


-54.878210 


20.45 ±0 


01 


1 


,87 ±0.06 


0.68 ± 0.02 




14086 


357.829155 


-54.897455 


20.48 ±0 


01 


1 


,68 ± 0.04 


0.59 ± 0.01 




18418 


357.768333 


-54.801860 


20.49 ±0 


.01 


1 


,65 ±0.04 


0.59 ±0.01 




13764 


357.905657 


-54.906287 


20.50 ±0 


01 


1 


,73 ±0.05 


0.65 ± 0.01 




10692 


357.819108 


-54.981725 


20.53 ±0 


02 


1 


,71 ± 0.06 


0.64 ± 0.02 




15516 


357.931958 


-54.867137 


20.56 ±0 


.01 


1 


,68 ±0.05 


0.61 ± 0.01 


0.3785 ± 0.0001 


19588 


357.961069 


-54.776725 


20.57 ±0 


.02 


1 


,64 ± 0.05 


0.64 ± 0.02 




15002 


357.897311 


-54.877856 


20.58 ±0 


01 


1 


,70 ±0.05 


0.67 ± 0.02 


0.3838 ± 0.0004 


14800 


357.901708 


-54.880859 


20.59 ±0 


.01 


1 


,73 ±0.05 


0.64 ±0.01 




15788 


357.914426 


-54.860804 


20.61 ±0 


.01 


1 


,75 ± 0.05 


0.64 ± 0.02 


0.3821 ± 0.0002 


15373 


357.965957 


-54.874282 


20.61 ±0 


01 


1 


,66 ±0.05 


0.64 ± 0.02 


0.3856 ± 0.0005 


13697 


357.905543 


-54.907479 


20.64 ±0 


.01 


1 


,81 ±0.06 


0.68 ± 0.02 


0.3782 ± 0.0005 


15187 


357.874035 


-54.873868 


20.65 ±0 


.01 


1 


,89 ± 0.06 


0.64 ± 0.02 


0.3868 ± 0.0006 


18026 


358.056024 


-54.810258 


20.66 ±0 


01 


1 


,74 ±0.05 


0.70 ± 0.02 




14378 


357.875627 


-54.892067 


20.71 ±0 


02 


1 


,66 ± 0.05 


0.56 ± 0.02 




14844 


357.899766 


-54.880469 


20.72 ±0 


.04 


1 


,75 ±0.31 


0.72 ± 0.09 




17455 


357.997121 


-54.822676 


20.72 ±0 


.04 


1 


,54 ±0.40 


0.88 ±0.12 




17729 


357.875358 


-54.816995 


20.74 ±0 


01 


1 


,78 ± 0.06 


0.64 ± 0.02 




15068 


358.094352 


-54.876409 


20.75 ± 


.02 


1 


,64 ±0.07 


0.68 ± 0.02 




15994 


357.763211 


-54.855887 


on qo _l_ n 


no 


1 


,80 ±0.06 


0.64 ± 0.02 




12892 


358.014272 


-54.926461 


on ofi _l_ n 
20.86 ± 


02 


1 


,65 ± 0.07 


0.68 ± 0.02 




15697 


357.897819 


-54.862968 


on Qa _l_ n 


no 


1 


73 ± 0.06 


0.55 ± 0.02 


0.3789 ± 0.0003 


12589 


357.785585 


-54.933352 


20.87 ± 


.03 


1 


66 ± 0.10 


0.63 ± 0.03 




11976 


357.893073 


-54.951394 


20.88 ±0 


02 


1 


62 ± 0.06 


0.61 ± 0.02 




14664 


357.901167 


-54.885034 


20.89 ±0 


.02 


1 


,84 ± 0.08 


0.68 ± 0.02 




13901 


357.824131 


-54.902094 


20.90 ±0 


.02 


1 


,69 ±0.06 


0.64 ± 0.02 




14595 


357.914211 


-54.886290 


20.93 ±0 


.03 


1 


,80 ±0.13 


0.66 ± 0.03 




12863 


357.874307 


-54.926847 


20.95 ±0 


.02 


1 


,66 ±0.07 


0.62 ± 0.02 




15156 


357.891357 


-54.874854 


20.95 ±0 


.03 


1 


,59 ±0.09 


0.56 ±0.03 




14825 


357.922988 


-54.880749 


20.95 ±0 


02 


1 


,63 ±0.07 


0.64 ± 0.02 




12650 


357.958048 


-54.931820 


20.96 ±0 


02 


1 


,77 ±0.08 


0.64 ± 0.02 




14407 


357.917097 


-54.891402 


20.97 ±0 


.02 


1 


,75 ±0.07 


0.63 ± 0.02 


0.3799 ± 0.0004 


14939 


357.872276 


-54.878410 


20.97 ±0 


02 


1 


,75 ±0.07 


0.62 ± 0.02 




14944 


357.913836 


-54.877816 


20.98 ±0 


.03 


2 


17 ±0.23 


0.64 ± 0.05 




Other Spectroscopic Cluster Members 














14271 


357.899268 


-54.896523 


19.13 ±0 


.01 


1 


,90 ±0.02 


0.71 ± 0.01 


0.3814 ±0.0002 


15403 


357.908847 


-54.870362 


19.52 ±0 


01 


1 


,58 ±0.02 


0.67 ± 0.01 


0.3860 ± 0.0003 


15827 


357.957606 


-54.860595 


19.64 ±0 


01 


1 


11 ±0.02 


0.45 ± 0.01 


0.3900 ± 0.0004 


15400 


357.866255 


-54.870713 


19.69 ±0 


.01 


1 


,62 ±0.03 


0.65 ±0.01 


0.3768 ± 0.0003 
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Table 3 — Continued 



Object ID a 


RA b 


Dec b 


i(MAG_AUT0) a 


g — r a 


r _ ja 


rcdshift z c 


14466 


357.909595 


-54.890780 


20.47 ±0.02 


1.57 ±0.05 


0.60 ±0.02 


0.3827 ± 0.0002 


14492 


357.914762 


-54.888447 


20.88 ±0.02 


1.58 ±0.06 


0.64 ± 0.02 


0.3899 ± 0.0004 


13372 


357.917399 


-54.914403 


20.98 ±0.02 


1.60 ±0.06 


0.65 ± 0.02 


0.3803 ± 0.0003 


14505 


357.870029 


-54.888283 


21.27 ±0.02 


1.63 ±0.08 


0.66 ±0.03 


0.3860 ± 0.0003 



a Objcct ID numbers are from the SExtractor catalog obtained using the i-band image for object de- 
tection. The objects are ordered from bright to faint by i-band MAGJVUTO, starting with the BCG. g — r 
and r — i colors are computed from 3"-diameter aperture magnitudes. The errors are simply statistical 
errors reported by SExtractor. Not included are photometric calibration errors estimated to be 0.03-0.05 
mag per filter. 

b RA and Dec are epoch J2000.0 and are given in degrees. 
c Redshifts measured from CMOS spectroscopy f 8 13.31 . 

d Galaxies, with i < 21, determined to be cluster members using maxBCG color selection criteria. 
Members are also limited to be within a radius t|qq = 1.51 h^ 1 Mpc (= 6.88') of the BCG. See ^13.21 for 
details. 

e Additional galaxies determined to be cluster members via GMOS spectroscopic rcdshifts f i|3,3l l, but 
which did not meet the maxBCG color selection criteria. 



-45 - 



Table 4. Other Galaxies 



Object ID b 


RA C 


Dcc c 


i(MAG_AUT0) b 




9 


_ r b 


r - i° 




redshift 


z d 


14193 


357.895093 


-54.901998 


17.99 ±0.00 


1 


.59 


±0 


01 


0.58 ± 0.00 


0. 


2970 ± 0. 


0003 


15313 


357.893518 


-54.875789 


18.75 ±0.00 


1 


11 


±0 


01 


0.47 ±0.01 





2486 ± 0, 


0002 


16682 


357.903652 


-54.840904 


19.97 ±0.01 


1 


.40 


±0 


03 


0.57 ±0.01 


0. 


3259 ± 0. 


0002 


13520 


357.887606 


-54.911328 


20.10 ±0.01 





.51 


±0 


01 


0.27 ±0.01 


0. 


0649 ± 0. 


0001 


19352 


357.902529 


-54.852090 


20.18 ±0.01 


1 


01 


±0 


02 


0.25 ± 0.01 





4214 ±0. 


0003 


15509 


357.941448 


-54.869154 


20.29 ±0.01 


1 


.57 


±0 


.05 


0.62 ± 0.02 





4178 ±0. 


0005 


16409 


357.890133 


-54.846245 


20.30 ±0.01 





.89 


±0 


.02 


0.31 ± 0.01 





3251 ±0. 


0002 


16570 


357.911876 


-54.843091 


20.42 ±0.01 





.77 


±0 


.03 


0.58 ±0.02 





1277 ±0. 


0002 


13423 


357.960803 


-54.913826 


20.63 ±0.02 


1 


01 


±0 


.03 


0.68 ± 0.02 





2524 ±0. 


0002 


13620 


357.889293 


-54.909472 


20.86 ±0.02 


1 


.66 


±0 


.08 


0.90 ± 0.02 





5354 ± 


0004 


19257 


357.902437 


-54.851578 


21.45 ±0.03 





.98 


±0 


.03 


-0.06 ±0.02 





2970 ± 0. 


0004 


16562 


357.948746 


-54.841891 


21.58 ±0.03 


1 


.90 


±0 


.12 


0.61 ± 0.03 





3595 ± 0. 


0002 



a Galaxies determined to be non-cluster members based on CMOS spectroscopic rcdshifts f i|3.3p . 

b Object ID numbers are from the SExtractor catalog obtained using the i-band image for object detection. 
The objects are ordered from bright to faint by i-band MAG_AUT0. g — r and r — i colors are computed from 
3"-diamctcr aperture magnitudes. The errors are simply statistical errors reported by SExtractor. Not 
included are photometric calibration errors estimated to be 0.03-0.05 mag per filter. 

C RA and Dec are epoch J2000.0 and are given in degrees. 

d Redshifts measured from CMOS spectroscopy ( £13.311 . 
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Table 5. NFW Fit Results 



Filter 



mi 



TOj'X (Qxx + Qyy^)r. 



z cr u d M 2OO (1O 14 M ) C C200 c X 2 /dof f 



tangential shear 














9 


lotto 


ZZ.O 


z4.U 


lo. / O 


U.Do 


U. 1 


r 


7013 


22.0 


24.0 


12.0 


0.70 


3.9: 


i 


3296 


22.0 


23.5 


12.0 


0.71 


5.9; 


z 


2300 


20.5 


22.5 


12.0 


0.62 


11.0 


i+r 


7995 








0.70 


1.2; 


i+r+z 


8996 








0.70 


5.0: 


i+r+z+g 


9424 








0.70 


1.3; 


j+r+ 2 +SL(s) g 


8996 








0.70 


4.8; 


i+r+2+SL(sg) g 


8996 








0.70 


1.9; 


WL+SL+cTc+Afaoo 11 


8996 








0.70 


5.1: 


B-mode shear 














9 


1883 


22.5 


24.0 


18.75 


0.68 


1.6: 


r 


7013 


22.0 


24.0 


12.0 


0.70 


o.r 


i 


3296 


22.0 


23.5 


12.0 


0.71 


0.1: 


z 


2300 


20.5 


22.5 


12.0 


0.62 


5.5^ 


i+r 


7995 








0.70 


o.r 


i+r+z 


8996 








0.70 


o.r 


i+r+z+g 


9424 








0.70 


o.r 



+0.4 

0.1 
+2.9 
2.1 
+5.3 
-3.8 
+ 11.9 
-7.1 



+2.8 
2.1 
+2.9 
2.3 
+2.8 
2.2 
+2.8 
2.2 
+2.9 
-2.2 
+ 1.3 
-1.3 



+3.1 
-1.5 
+0.1 
0.1 
+0.4 
0.1 



•5.2 



+ 0.1 
-0.1 
+0.1 
0.1 
+0.6 
0.1 



> 45 





83 


0. 


68 


6 5+ 5 ' 3 
D -°-3.0 


1. 


54 





059 


7 +13.1 
°- ' -2.6 


1 


55 


0. 


055 


i - 8 -1.8 





89 





60 



e.i±S:8 


1 


.58 





,048 




1 


.48 





.077 




1 


.50 





,069 


6.2+^7 


1 


.48 





,077 


5 5+ 2:7 


1 


.48 





,077 




1 


.48 





,077 



B 5+ 10 1 
°-°-5.4 


1.04 


0.41 


> 63 


1.19 


0.25 


> 


0.91 


0.58 


3+ 1 ' 2 


0.61 


0.91 



> 51 


1.10 


0.34 


> 27 


0.80 


0.71 


> 


0.77 


0.75 



a Numbcr of galaxies used in the weak lcnsing analysis in each filter. For the multi-filter samples N is the number of unique 
galaxies. 

b SExtractor MAG_AUT0 magnitude limits used to define the galaxy sample. 
c Minimum Ellipto size Q X x + Qyy used to define the galaxy sample. 

d The source redshift at which 1/S cr it is the same as the effective value computed by integration over the source galaxy 
redshift distribution, as described in 35.31 

c Best-fit NFW profile parameters: mass M200 an d concentration c 2 oo- Errors are 1-paramcter, la values, as determined 
by where Ax 2 = 1. The uncertainties on Af 2 oo arc rounded off to the nearest 0.1 X 1O 14 M0. Note that for most of the 
cases (primarily B-mode fits) where there is no significant mass detection, we provide only a lcr lower limit on c 2 oo, which is 
otherwise not constrained on the high side even at lcr, up to the upper bound value c 2 oo = 10 4 that we have checked. Joint 
2-parameter error contours for select samples are shown in Fig. 1131 

f X 2 P er degree of freedom (dof) relative to a null hypothesis of zero shear. (This is not the x 2 /dof of the NFW fit, which 
is very close to one in all cases.) P is the probability of exceeding the observed x 2 /dof. The number of degrees of freedom 
for this x 2 t es t is always 20, i.e., the number of radial bins plotted in Figurcs llOH2l 

g Fit results derived from combined weak plus strong lensing ("SL") constraints, "(s)" denotes the case where we estimated 
the dark matter mass within the Einstein radius by subtracting off just a stellar mass contribution, while "(sg)" is the case 
where we also subtracted off an estimated gas mass contribution. Sec *J67TJ for details. 

h Fit results derived from combined weak lensing (i + r + z), strong lensing (SL(sg)), cluster velocity dispersion (a c ), and 
cluster richness (N200) constraints. See 36.21 for details. 



Table 6. Source Galaxy Star Formation Rates' 



Knot 


fiy)[Oll] (erg s 1 cm 2 ) 


/Mi (erg s 1 cm 2 ) 


f(v)g (erg s 1 cm 2 Hz 1 ) 


SFR ((/ ie „ s = 39) M yr" 1 ) 


SFR ((/ iens = 141) M yr" 1 ) 


Al 


1.06 ±0.04 x 10~ 15 


1.71 ±0.06 x 10- 16 


1.36 ± 0.02 x 10~ 28 


3.9 ± 1.1 


1.1 ±0.4 


A2 


0.84 ±0.04 x 10" 15 


1.43 ±0.06 x 10~ 16 


1.43 ± 0.02 x 10~ 28 


3.1 ±0.9 


0.85 ± 0.4 


A3 


2.09 ±0.10 x 10~ 15 


1.28 ±0.06 x 10~ ie 


0.51 ± 0.01 x 10~ 28 


7.7 ±2.2 


2.1 ±0.4 


Al 


1.02 ±0.02 x 10~ 15 


2.83 ±0.06 x 10- ie 


2.33 ± 0.02 x 10~ 28 


3.7 ± 1.1 


1.0 ±0.4 



a See J7]for the definitions of the various fluxes f(u). Fluxes quoted are measured values. fi ena is the lens magnification. 



